The magnetoelectric coupling in the multiferroic compound L1CU2O2 
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We investigate the possible types of coupling between ferroelectricity and magnetism for the 
zig-zag spin chain multiferroic LiCu^O^ compound. We construct a multi-order parameter phe- 
nomenological model for the material based on a group theoretical analysis. From our calculation 
we conclude that a coupling involving the inter-chain magnetic structure and ferroelectricity is nec- 
essary to understand the experimental results of Park et.al. [J]. Our proposed model is able to 
account for the electric polarization flop in the presence of an externally applied magnetic field. 
Furthermore, based on our theoretical model we can make specific selection rule predictions about 
electromagnon excitations present in the L1CU2O2 system. We also predict that the electromagnon 
peaks measured in an ac-conductivity measurement are field dependent. 
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I. INTRODUCTION 

Cross coupling of magnetism and ferroelectricity in 
material systems is an intriguing phenomenon. Presently 
there are chemical compounds in which both magnetism 
and ferroelectricity can exist simultaneously. These novel 
systems are called multiferroics @, 0, 0, H,IE 0j- 

Multiferroic materials have the unique feature that 
their magnetic properties can be controlled by an elec- 
tric field and the electrical properties by a magnetic field. 
Such a control allows for the possibility to fabricate mul- 
tifunctional devices with great potential for technological 
applications in areas such as memory, sensor, and spin- 
tronic devices @, [t| [l(| [H| • The present interest in multi- 
ferroics arises from the controlled tunability of these sys- 
tems. Examples of multiferroic systems include the rare 
earth perovskites RMnOs. RMn 2 0§ (R: rare earths - 
Gd, Dy, Tb...) [lUmEljll, the delafossite CuFe0 2 
l|, the spinel CoCr 2 4 [3], MnW0 4 [H, Ni 3 V 2 O s 
18|, the hexagonal ferrite A 2 Zn 2 Fe\ 2 22 (A = Ba,Sr) 
19 1, the zig-zag spin chain compound LiCu 2 2 [H and 



The mixed valent zig-zag spin chain cuprate compound 
LiCu 2 2 has been shown to be a multiferroic material by 
Park et. al. [l[ • The experimental data on the compound 
exhibits the canonical multiferroic behavior of a 7r/2- flip 
in the direction of its electrical polarization for a mag- 
netic field perpendicular to the axes of polarization. The 
data also reveals that the compound exhibits a unique 
behavior of generating an electric polarization parallel 
to the applied magnetic field direction. Due to sample 
issues, the second experimentally observed feature has 
to be investigated further before being confirmed exper- 



imentally [20(. Recently, a theoretical attempt has been 
made to explain these behaviors for the LiCu 2 2 com- 
pound using a stoichiometry argument [U . In this 
article, we will not provide any physical explanation for 
the occurence of this effect within our proposed model 
and will focus on explaining the polarization flip only 
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In the microscopic theories |23|, |24|, |25( the polarization 
flip has hardly been studied so far. An explanation of the 
7r/2- polarization flip in the conventional continuum the- 
ory [261 ] requires an introduction of numerous anisotropic 
magnetic terms which is rather unnatural. In this paper 
we take a different approach using both a group theo- 
retical analysis involving the symmetry of the magnetic 
structure [271 ] and a phenomenological multi-order pa- 
rameter magnetization model to study the experimental 
data [l[ . The multi-order parameter model is able to ac- 
count for the electric polarization flip through 7r/2- in 
the presence of an applied external magnetic field. Wc 
show that a coupling involving the inter-chain magnetic 
structure and ferroelectricity can be used to understand 
the experimental results [l[ and in particular the polar- 
ization flip. The theory does not include any anisotropic 
terms. Furthermore, the presence of inter-chain coupling 
is supported by experimental evidence found in the Ra- 
man scattering experiment (28j . 

The new inter-chain coupling term cannot be expressed 
in the familiar P ■ [(M • V)M - (V • M)M] form in the 
continuum limit. A hidden assumption behind the con- 
tinuum coupling model [2rJ| is that the magnetic order 
is described by a single order parameter. When the lat- 
tice structure is complicated this assumption is not valid 
and the general form of the magnetoelectric coupling can 
be different. The additional inter-chain magnetoelectric 
coupling derived for the LiCu 2 2 compound is a reflec- 
tion of the physics that there are novel magnetoelectric 
coupling terms in a multi-order magnetic structure. Be- 
sides accounting for the observed experimental features 
we also discuss the selection rules associated with the low 
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energy electromagnon excitations. 

This paper is arranged in six sections. In section Til Al 
we describe the crystal and magnetic structure of the 
L1CU2O2 compound under study. In section III Bl we 
discuss the experimental evidence which indicates that 
inter-chain coupling has an important role to play. In 
section IIIII we elucidate the magnetic symmetry opera- 
tions of the lattice and use group theory to construct the 
possible intra- and inter- chain couplings present in this 
compound. In section Hvl we propose the phenomenolog- 
ical model for the L1CU2O2 system. In section IVl we first 
state and then analyze the theoretical model which ex- 
plains the experimental data of Park et.al.l]. We then 
derive the selection rule associated with the hybrid exci- 
tations of phonon and magnon termed electromagnons. 
Finally, in section I VII we provide a summary and state 
the main conclusions of our paper. 



II. THE LiCu 2 2 SYSTEM 

A. Crystal and magnetic structure 

The mixed valent cuprate compound L1CU2O2 is a fer- 
roelectric material. It crystallizes in an orthorhombic 
structure with the space group Pnma (# 62) (29l[30l. [3l| . 
The crystal structure of the compound can be visualized 
as follows. Consider two linear Cu 2+ chains propagat- 
ing along the crystallographic 6-axis. The two chains 
are displaced b/2 with respect to each other and they 
form a zig-zag triangular ladder like structure as shown 
in Fig. [T] The ladders are separated from each other by 
the non-magnetic Li + ions (distributed along the crystal- 
lographic a-axis) and by the layers of nonmagnetic Cu + 
ions (distributed along the crystallographic c-axis). The 
magnetic behavior of the system is provided by the Cu 2+ 
ions which carry a spin-1/2. The letters a, 6, and c rep- 
resent the lattice constants of the simple orthorhombic 
crystal structure. 




FIG. 1: Spin chain: The multiferroic LiCu-zO-z compound has 
pairs of Cu 2+ chains running parrallel to each other. Each 
chain is separated from the other by b/2 and forms the zig- 
zag triangular ladder structure shown in the figure. The black 
dots in the figure represent the magnetic Cu 2+ ions. These 
ions carry a spin-1/2 indicated by the arrow in the figure. The 
ladder is directed along the crystallographic 6-axis. 

The magnetic structure of LiCu 2 02 has been deter- 
mined from the neutron scattering experiments [32| ■ The 
magnetic modulation vector obtained from these experi- 
ments is Q = (0.5^, 0) where f = 0.174 is the spiral 
modulation along the chain direction. The letters a, b, 



and c represent the lattice constants of the simple or- 
thorhombic crystal structure as mentioned earlier. Along 
the a-axis there is antiferromagnetic order and along the 
6-axis there is a spiral order. Succesive spins on each rung 
are almost parallel to each other with each spin being ro- 
tated relative to each other by an angle a = 2-7r£. Within 
each double chain any nearest-neighbor spins from op- 
posite legs are almost antiparallel and form an angle 
a/2 = 7r£. Neutron diffraction experiments [32j conclude 
that rotating spins lie in the be plane [33| . 



B. Experimental evidence of inter-chain coupling 

The present theoretical formulation of spiral ferro- 
electrics do not consider the role of inter-chain coupling. 
For the L1CU2O2 compound, a spiral ferroelectric, there 
is experimental evidence from the Raman scattering ex- 
periment of Choi et. al. [28[ which suggests the presence 
of inter-chain interactions. Referring to Fig. 1 of their 
paper one can clearly observe that similar peaks are ob- 
served for both the parallel (xx) and crossed polariza- 
tions (xy). This indicates that inter-chain effects are 
important and should be considered in the theoretical 
formulation of LiCv,2C>2- Similar behavior has also been 
observed in the Na-analog compound NaCu 2 02 [34ll35l|. 
In this article we provide a theoretical justification using 
group theory for this important experimental observa- 
tion. We also construct a phenomenological theory con- 
sidering the effects of inter-chain coupling in multiferroic 
compounds. 



III. SYMMETRY OPERATIONS OF THE 
LATTICE AND MAGNETIC STRUCTURE 

The goal of this section is to find a physically appro- 
priate form for the magnetoelectric coupling in L1CU2O2 
subject to the symmetry constraints of the magnetic 
structure. Physically such a term should conserve trans- 
lational symmetry, space inversion symmetry, and time 
reversal symmetry. Phenomenological models (2(| are 
based on these considerations only. These models as- 
sume that a single unit cell of the lattice is a point in the 
system without any structure. An appropriate realistic 
model in these cases should include all the inequivalent 
atoms (considering both the type and location of the ion 
in space) in a single magnetic unit cell. This allows for 
the possibility to include more degrees of freedom. For 
example, in the LiCu202 case there are four spatially in- 
equivalent magnetic ions in the magnetic unit cell that 
we should consider. They are denoted as Si, S2, S3, and 
S4 respectively (see Fig. [2]). 

As the number of possible forms of the magnetoelec- 
tric coupling increases rapidly with the degrees of free- 
dom such a detailed model seems too complicated to be 
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FIG. 2: The dashed and the bold line indicates the bond be- 
tween the Cu 2+ ions forming the zigzag chain, (a) Magnetic 
unit cell. The four black dots represent Cu 2+ ions, (b) Side 
view of the same magnetic unit cell. The center of inver- 
sion is situated halfway between the two Cu 2+ ions. These 
ladders are separated from each other by the non-magnetic 
Li + ions (distributed along the crystallographic a-axis) and 
by the layers of nonmagnetic Cu + ions (distributed along the 
crystallographic c-axis) . 



useful. However there are symmetry constraints which 
help to simplify the situation. These constraints should 
be obeyed by any term in the Hamiltonian including 
the magnetoelectric coupling term. The terms should 
be invariant under the lattice symmetry operations and 
time reversal symmetry. The lowest order magnetoelec- 
tric coupling which preserves time reversal symmetry is 
a trilinear term. It is linear in the electric polarization 
and bilinear in the magnetic order parameters. To sim- 
plify the discussion, we consider the coupling terms that 



involve the uniform spontaneous electric polarization, P, 
and not those including the modulations of the polariza- 
tion. The overall polarization has been experimentally 
measured and confirmed. The general magnetoelectric 
coupling is 

H ME = A^P«^£aW£-£0Sf(£)S7(i?')(l) 

a/37 ij RR> 

where a, ft, 7 = x, y, z; i, j = 1, 2, 3, 4, and 
R denotes the position of the unit cell. The magnetic 
structure in the compound can be written as S?(R) = 
S" (Q) exp(iQ ■ R) + c.c. Substituting this expression into 
the above Hamiltonian we have 



H 



ME 



a/37 y 



c.c) (2) 



where Cij ;ct ^ 7 ((5) is the Fourier transform of CV, ;Q( g 7 (-R — 
R>). 

Not every term in the above magnetoelectric coupling 
conserves all the symmetries. Only certain linear combi- 
nations do. Specifically we are interested in the lattice 
symmetry operations that conserve the magnetic struc- 
ture with the magnetic propagation vector Q. From 
group representation theory we know that these oper- 
ations constitute a subgroup of the full symmetry group 
of the lattice for which we have one (and only one) two- 
dimensional representation, E, 2h,m c and m c 2h, all being 
two-by-two matrices . The lattice symmetry operations 
which preserve the magnetic structure and the magnetic 
propagation vector Q are listed in Table Q] 



TABLE I: Lattice symmetry operations which preserve the 
magnetic structure and in turn the magnetic propagation 
vector- Q. The identity operation is represented by E. A 
two-fold rotation about the crystallographic 6-axis is repre- 
sented by 2^. A reflection about the c-axis is denoted by 
m c and finally a combination of the the rotation-reflection by 
m c 2 b . 

Er = (x, y, z) 

2 b r=(x,y + l/2,z) 

m c r = {x-l/2,y,z-l/2) 

m c 2 b r = (x- 1/2, y + 1/2, z - 1/2) 



The two dimensional representations are given by 

,— iqb/2 



E 
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;2 6 = ( 
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1 



: m r 2h 



e 






-iqb/2 



_ e -iqb/2 





(3) 



(4) 



For the above representation we can find groups of sym- 
metry adapted variables which are transformed accord- 
ing to the representation under the symmetry operations. 
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These variables can be constructed as follows: a symme- 
try adapted variable is a column vector (two-dimensional 
in this case) whose elements are a linear combination of 
the original magnetic variables - Sf(Q)'s. The variables 
transform under any operation of the subgroup in the 
same way they transform if left-multiplied by the corre- 
sponding matrix for the same operation. These variables 
can be computed for the LiCu 2 2 case after deducing 
the transformation tables (see Table |TT] and Table IIIip of 
all the members of the symmetry subgroup. 

TABLE II: Symmetry operation: 2b, two- fold rotation about 
the crystallographic 6-axis. In LiCu2 02 there are four spa- 
tially inequivalent magnetic ions in the magnetic unit cell. 
They are denoted as Si, S2, S3, and 54 respectively (see 
Fig. [2}. The subscripts a,b, and c denote the components 
along those crystallographic axis. 

Si S2 S3 S4 



Sig — —Si a e~ lqb ^ 2 Szg = —S4 a e~ lqb / 2 S Sa = —Si a S 4a = — 52a 
S lb = — 5*3be lqb / 2 S 2 b = Site lqb ^" S sb = Sib S 4b = 8 2b 
Si c = —S3 C e lqb / 2 S 2 c = —S4, c e lqb / 2 Sg c = —Sic S4C — —S2c 



J 



Using the properties of the symmetry adapted variables 
I 



Hereafter we suppress the argument of S. We will con- 
sider S to be S(Q) and S + to be S(—Q). Having found 
the symmetry adapted variables, we put the magneto- 
electric coupling in the following general form 

Hme = A P a S f^ M H^ S ^ + c - c - ( 9 ) 



As a result six groups of symmetry adapted variables 
are found, that is, twelve elements in total. The number 
of symmetry adapted vectors can also be deduced if we 
notice the fact that all of their elements should make a 
new basis for the magnetic structure which requires the 
total number of elements to be twelve, the same as the 
original S"'s, and the total number of vectors to be six. 
We divide them into two sets as listed below. The groups 
of symmetry adapted variables are 



TABLE III: Symmetry operation: m c , reflection about the 
crystallographic c-axis. In LiCu2 02 there are four spatially 
inequivalent magnetic ions in the magnetic unit cell. They 
are denoted as Si, S 2 , S3, and S4 respectively (see Fig. 
The subscripts a, b, and c denote the components along those 
crystallographic axis. 



Si 


s 2 


S 3 


s 4 


Sla = —S2a 


S2a = Sla 


S3a = S4a 


S4a = — S3a 


S lb = — S26 


S26 — Slb 


S36 = S46 


S45 = — S36 


Sl c = S2c 


S2c = —Sic 


S3 C = — S4c 


S4 C = — S3c 



(5) 



(6) 

I 

we have 



= sf\a z M l3 . l3l cj z )S ( j ) (7) 
= Sf\a v M ij . M )S^. (8) 

I 

where i, j — 1, 2 and M is an arbitrary two-by- 
twocoupling matrix. Using the properties of the sym- 
metry adapted variables we have 



Si, + e-^ 2 S 3a 
S 2a - e-"* b / 2 S ia 



S lb - e-^ b / 2 S 3b 
S 2b + e-^ b ' 2 S ib 



S 2c 



e -iqb/2 Si 



S 2a + e^ b / 2 S 4a 
-S la + e-^Sia 



526 



-s 



lb 



-iqb/2 g 
-iqb/2g : 



lb 



:S/> 



Slc + e- l ^ 2 S 3c 
-S 2c + e~^ b / 2 S ic 
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Sf ] M^S^ ^ Sf\m 2b y M ij;f3l m 2b S^ = (a z M ij;01 a z )S^ j) (10) 
Sf^M^S^ ^ S ( ^(m mc yM ir ^m mc S^ = sj? f (a y M ir ^a y )S^ (11) 



Now we can apply the stated constraints to find all the 
possible forms of the magnetoelectric coupling. 

1. Time reversal symmetry: The expression above 
with its tri-linearity automatically includes time re- 
versal symmetry. 

2. Lattice symmetry: We focus on the subgroup of 
the lattice symmetry operations. All three compo- 
nents of the electric polarization behave differently 
under the symmetry operations. We discuss each 
case separately. In the following we explicitly work 
out the case for P c and simply list the result for 
the other two components. The symmetry prop- 
erties of P c under the lattice symmetry operations 
arc P c ^ -P c , P c ^ -P c , P c -P c . From the 
overall invariance of the trilinear coupling we re- 
quire that Mij./3j anti-commute with both a y and 
a z . From this we can infer that My should be pro- 
portional to a x to preserve the invariance under 
symmetry operations. A similar procedure can be 
applied to find the appropriate My's for P a and 
Pb- The result can be summarized as follows: for 
P a , My (x a y ; for P b , M {j cx 1; for P c , My oc a x . 
Later we will prove that the constant of proportion- 
ality can be either real or purely imaginary based 
on space inversion symmetry arguments. 

3. Inversion symmetry: Inversion operation is not a 
I 



if (3, 7 = a, c or [3 = 7 = 6, and 

P c S^a x S^ ±{-P c )(-Sf T o- z a,. 



member of the subgroup of the symmetry opera- 
tions that conserve the magnetic propagation vec- 
tor Q. Therefore it cannot be represented by a 
matrix acting on the symmetry adapted variables. 
The mapping for space inversion operation is 

S?(q) - SS(-q), S?{q) - S2(-q) (12) 
S3 (?) - Sn-Q), S2(q) - S?(-q) (13) 

where a denotes the spin component. With this 
mapping we obtain the expression for the inversion 
operation in terms of the symmetry adapted vari- 
ables as 



where a = a,c and i = 1, 2 are the indices for the 
symmetry adapted variables. Below we explicitly 
work out the case for P c and list the result for the 
other two components of electric polarization. We 
have 



■ z S^*) = P c (S^a x S^y, (16) 



. zS U>) = -P c (sMa x spr, (17) 



if one and only one of (3 and 7 is b. The elec- not violate any of the symmetries of the system 

trie polarization and the total Hamiltonian must 
be real making the proportionality constant in the 
first case real and in the second case purely imagi- 
nary. The complete result is summarized below. It 
shows that there are only six generic terms that do 



P a (S^a y S^ ] +c.c) (18) 
P b (S^S^+c.c) (19) 
P C (S { ; )] 'a x S^ + c.c) (20) 
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if (3, 7 = a, c or f3 — 7 = 6, and 



adapted variable as 



iP b (S^S^ - 



<J) 



— c.c) 
c.c) 



c.c 



if one and only one of (3 and 7 is 6. 



(21) 
(22) 
(23) 



To summarize our work up to this point we have ex- 
hausted all possible forms of magnetoelectric coupling 
that are invariant under (1) time reversal, (2) space in- 
version, and (3) all the lattice symmetry operations that 
conserve the magnetic propagation vector by using the 
symmetry adapted variables. 



IV. PHENOMENOLOGICAL MAGNETIZATION 
MODEL 

A simplified expression for the phenomenological 
model can be obtained if we observe first that ferroelec- 
tricity coexists with the non-collinear magnetic structure. 
This suggests that terms with (3 = 7 can be excluded. 
Second, the polarization along the 6-axis is not observed. 
Therefore, the couplings with P b need not be considered. 
Third, the magnetic moments on the four Cu 2+ atoms, 
assumed to be independent variables in our theoretical 
formulation, form two zigzag chains [1^, [3(J, HH extended 
in the 6-direction on each of which a spin density wave 
with a propagation vector q along the fe-axis exists. To 
be specific, in Fig. [1] they propose that S2 and S4 are on 
one zigzag chain and Si and S3 in the adjacent unit cell 
are on another zigzag chain. This observation drastically 
simplifies the expression of magnetoelectric coupling as 



S a = e -iqb/2 S a 

S% = -e-^^St 



(24) 
(25) 



These relations reduce the symmetry adapted variables 
to their final form as displayed in Eq. 1281 Considering 
the fact that S 2a = -e- lqb / 2 S ia and S u = e' iqb / 2 S 3b 
the symmetry adapted variables S^ 1 ) and iS^ 2 ) become 



5(2) 



2e~ lqb / 2 S 3a 
-2e~^ 2 S Aa 



-2e~ lqb / 2 S 4b 
-2 e - iqb / 2 S 3b 



2e~^ 2 S 3c 
2e- ic > b l 2 S ic 



(26) 



(27) 



All the magnetoelectric couplings terms can be enumer- 
ated in terms of the symmetry adapted variables by using 
the above two equations. After the insertion of Eqs. [2H 
and [53] three of the six symmetry adapted variables re- 
duce to zero. We can instead take the reduced symmetry 



S = 



c(l) 
o(2) 



,(2) 
6 

.(1) 



>t 3 



(28) 



where the superfix denotes the chain number to which the 
spin belongs and the argument of the <S"s are suppressed. 
However one should bear in mind that all the «S"s are 
Fourier components and in general complex numbers. 



A. Intra-chain interaction terms 

First let us enumerate all the intra-chain couplings that 
couple to P c and P a . Since P b has not been observed 
under any conditions we ignore it in our theoretical con- 
sideration. These couplings are 



£P c (SpS«- c.c. ±(1)^(2)) 



and 



5«+c.c.T(l)-(2)) 



(29) 



(30) 



where the upper sign in the ± (=f) is for the case when 
the subscript is c. 

Our calculations show that an existing phenomeno- 
logical theory [26] for a spiral ferroelectric material in 
which only intra-chain interactions are included does not 
suitably explain the present system of interest. We de- 
mostrate through our calculations that some of the unex- 
plained experimental facts may be ascribed to inter-chain 
interactions generated as a result of considering the de- 
tails of the magnetic unit cell. We also show that not all 
the intra-chain and the inter-chain terms have the usual 
P ■ \{M ■ V)M - (V • M)M] phenomenological form. 

In the previous theory it is assumed that the electric 
polarization along the c-axis can be explained by its cou- 
pling with a cycloidal spin density wave on the 6c-plane. 
Referring to Eq. [201 one finds that P c must be coupled 
to S , 4f,S'| c + S 3b S 3c + c.c As we know from above that 
S4 and S 3 are the order parameters of different zigzag 
chains, we conclude that in this theory, the ferroelec- 
tric order P c couples to bilinear products inside a zigzag 
chain. Furthermore, when the magnetic field is applied 
along the 6-axis, P a can also be induced, and it is ex- 
plained by its coupling to a cycloidal spin density wave 
in the a6-plane. Using Eq. [18] we know that P a can only 
couple to S4 a Sl b + S 3a S^ b + c.c. which is a product inside 
a single zigzag chain. Following this line of reasoning one 
easily sees that at least some unexplained multiferroic 
properties of the system must come from magnetoelec- 
tric couplings that involves more than one chain, that is, 
an inter-chain coupling term. 

The physical picture of ferroelectricity induced by 
twisted magnetism has been widely recognized in mul- 
tiferroic theories [26[ . Through our calculation we point 
out that this physical picture is based on a single chain 
one-dimensional physics and has to modified to include 
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inter-chain interaction. This observation is based on sym- 
metry arguments and therefore promises to be indepen- 
dent of the specific form of interaction. 

In transforming the above expressions, which are the 
Fourier components, to an expression in the real space, 
one should notice that a imaginary unit i corresponds to 
an odd number of derivatives while the absence of an i 
corresponds to an even number of derivatives. Therefore, 
if we only consider the lowest order in the continuous 
limit the above expressions can be written as 

X c P c (S { b 1) (y)d y S^(y) - d v sj J 1 \y)Si%) 

±(1) -» (2)) 

(31) 

and 

P a (si 1) (y)S^ a (y) T (l)->(2)) (32) 

If the magnetic structure is a cycloidal spiral on the fee- 
plane we can easily see that P c ^ which is in accord 
with the phcnomcnological theory. The coupling in front 
of the P c component of polarization is dependent on 
the magnetic field along the a-axis and is sensitive to 
it. A discrepancy with the intra-chain phenomenlogi- 
cal theory occurs when we consider the observed ferro- 
electricity flop from the c-axis to the the a-axis. In the 
phenomenological theory it is explained as derived from 
a spin flop from the cycloid on the oc-plane to the ab- 
plane. However if we consider a magnetic structure such 
as S(x) = S (cos(qx) , sin(qx) , 0) we find that P a gives a 
modulated value and therefore averages to zero. Since the 
starting point from group theory is more general we be- 
lieve that the phenomenological expression for the mag- 
netoelectric coupling 

P ■ [(M ■ V)M - (V • M)M] (33) 

does not capture the entire physical picture. Specifically 
in this case the P C [{M -V)M -(V '• M)M] C has the correct 
form while the expression for P a [(M ■ V)M - (V • M)M] a 
is invalid atleast for the LiCuiOi case. The derivation of 
the phenomenological theory [2y] assumes certain general 



symmetries of the Hamiltonian to construct the magne- 
toelectric coupling for P a and P c . Based on such an argu- 
ment there is no reason for the coupling to be different. 

The starting point of the phenomenological theory [26| 
considers each unit cell as a single site with magnetization 
M. This assumption does not hold if the little group of 
the lattice symmetry operations has no one-dimensional 
representation. If the system can be effectively described 
by only the unit cell magnetization M, then all the little 
group operations should commute and there must exist 
at least one one-dimensional representation. In the case 
of multiferroics such as RMn0 3 0, [T3] and Ni 3 V 2 O s 
[l8| there are one-dimensional representations and the 
phenomenological theory provides a correct description 
[271 ] . For other multifcrroic systems such as i?Mn2C>5 
and L1CU2O2 only two-dimensional representations exist 
because of the antiferromagnetic structure along the a- 
axis. As a result the conventional phenomenological the- 
ory doesnot apply. From the above discussion we clearly 
see that the intrinsic structure of the unit cell is impor- 
tant and therefore must be taken into consideration. 



B. Inter chain interactions 

From the expressions for the intra-chain terms we can 
conclude that an explanation for the net magnetization 
along the a-axis cannot be obtained within the intra- 
chain theory. We therefore go beyond the scope of that 
theory and include the effects of inter-chain coupling. 
Since there are several terms involved we first focus on 
those that couple to P c . These expressions are 

P C (SW x SW) b (34) 
P c (X a S^S^ + X b S^S^ + \ C SPS1 2) ) (35) 

where the A's are arbitrary constants. Physically the 
first term could induce a polarization when the spins are 

aligned on the ac-plane and S'W-LS't 2 ). If the A's are 
close to each other the second term can be viewed as 
XP C (S± ■ S2). This term is minimized with a collinear 
spin arrangement. The terms that couple to P a are 



PaiPaS^dyS^ - (3 a dyS^S^ + PbS^dyS^ 



Same as above we can see that the first term induces 
a nonzero ferroelectricity with the alignment on the ac- 
plane and the collinearity between the two chains while 
the second term is minimized when SC 1 ) JlS^> . The above 
group theory analysis yields the possible intra-chain and 



P a (SW x dySW + X OyS^b (36) 
PbdyS^S^ + pjsPdySW - (3 c d y sPs^) (37) 

I 

inter-chain magnetoelectric terms. The coupling param- 
eters in the interactions cannot be determined without 
a proper understanding of the microscopic theory of the 
material compound and without improved experimental 
results. 
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V. EFFECTIVE HAMILTONIAN AND 
ELECTROMAGNON SELECTION RULES 

In this section we conjecture a simple effective Hamil- 
tonian which can explain the basic phenomena observed 
in the LiCuiOi experiments In particular we provide 
an explanation for the polarization flip in the presence of 
an external magnetic field. Furthermore, using the phe- 
nomenological Hamiltonian we also predict a selection 
rule to be obeyed by the hybrid excitations of phonon 
and magnon termed as electromagnons [36[ in the liter- 
ature. It is hoped that these electromagnon excitations 
can be experimentally detected in future experiments. 

A. Effective Hamiltonian 

The effective Hamiltonian that we propose includes 
only two terms. One term involves the intra-chain cou- 
pling which is responsible for the electric polarization 
along the c-axis. The other term involves the inter-chain 
coupling that is responsible for the electric polarization 
along the a-axis. The magnetoelectric coupling Hamilto- 
nian, H em , is given by 

H em = \ C P C (S^ X dy&V + X dyS^) a 

+\ a P a (SW x d y SW + £(2) x d y S(% (38) 

where the symbols have the same meaning as before. 
Both these terms satisfy the underlying lattice and mag- 
netic symmetry requirements. 

The proposed effective Hamiltonian explains the basic 
physics of the L1CU2O2 system. First, in the absence 
of any external magnetic field according to the neutron 
scattering experiments [32| the spins lie on the 6c-plane. 
We have for the ground state spin configuration, S^ and 
So , in the two chains 

§P = S(0,cos(gy),sin(gy)) (39) 
<?( 2 ) = 5(0, cos(qy + S(/)),sm{qy + 60) (40) 

where 5<fi is the phase difference between the spin con- 
figurations of the two chains. Based on Eq. [35] this spin 
configuration leads to an electric polarization along the 
c-axis 

P c cx (S (1) x dyS^ + S {2) x d y S {2) ) a = const. (41) 

Second, in the presence of an applied magnetic field along 
the 6-axis the spin should flip to the ac plane when the 
field is larger than a certain critical value. This spin-flop 
transition is expected in a very general magnetic model 
with relatively isotropic couplings and spiral magnetism. 
This is similar to the spin-flop transition in the antiferro- 
magnetic Heisenberg model in the presence of a Zeeman 
magnetic field. Therefore the expected ground state spin 
configuration in the two chains are 

= S(cos(qy),0,sm(qy)) (42) 
<?( 2 ) = S(cos(qy + 5(1), 0, sin(gy + 6<j))) (43) 



Based on Eq.[38]this spin arrangement leads to an electric 
polarization along the a-axis when the spin structure in 
both the spin chains are in phase, that is, 50 = 

P a cx (HP) x d y S {2) + S {2) x d y S {1) ) b = const. (44) 

The theoretical model developed above suggests that 
the polarization in the two crystallographic direction 
comes from two different types of coupling. The polariza- 
tion along the c-axis arises from the intra-chain coupling 
while the polarization along the a-axis stems from the 
inter-chain coupling. Furthermore, in formulating this 
model we do not include strong anisotropy in the spin 
model to explain the flop transition of the electric polar- 
ization in the presence of an applied external magnetic 
field. 

The proposed Hamiltonian explains the polarization 
flop by an applied magnetic field without introducing 
complicated anisotropic magnetic couplings [26]. The 
above model predicts explicit magnetic configurations 
corresponding to the polarization directions which can 
be verified in future experiments. 

B. Electromagnon selection rule 

Electromagnons are a combination of phonon and 
magnon excitations [36j . From the phenomenological 
model we can derive an explicit selection rule governing 
the electromagnons. The selection rule can be stated as 
follows: Electromagnetic waves that are polarized per- 
pendicular to the bulk electric polarization can be ab- 
sorbed. To be more specific in zero magnetic field only 
those electromagnetic waves that are polarized along the 
crystallographic a-axis can couple to the magnons. How- 
ever, with an applied magnetic field along the 6-axis only 
those waves polarized in the crystallographic c-axis can 
couple to the magnons. 

The selection rules can be obtained in the following 
manner. We first consider the case with the external 
magnetic field absent. The dynamics of the system can 
be derived by assuming a small deviation in the ground 
state properties 

P = P + u (45) 
§W = <?« + 5S W (46) 
S^ = S& + 5S& (47) 

where P , S X \ and S 2 ^ are the ground state values 
and the small deviations are indicated by u, 5S^ X \ and 
SS (2 \ Ener gy minimization requires that the first order 
terms vanish in the ground state. We therefore study 
the second order terms in the perturbation expansion to 
see how the dynamical degrees of freedom are coupled. 
Also the spin being a length-preserving vector we know 
that SSi-LSiQ. Defining rti = (0, cos(gy), sin(gy)), n,2 = 
(0, — sin(gy), cos(gy)), we have 5Si — m^r^ + m^a 
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(where a is the unit vector along the a-axis). We now 
consider the second order coupling that involves u c and 
SSi derived from Eq. M 

u c (5SW x d y S ( o ] - dySS^ x !$\ 

= S 2 u c (q(m^ n2 + m^a) x fi2 

— (9j,mpn2 - qm^ni + d y a) x ni) n 

= S 2 d y m 1 u c (48) 

In momentum space the above expression can be written 

as 



(49) 



From this we infer that the zero mode phonon u c {k = 0) 
does not couple to the magnons. Therefore, if we com- 
pute the optical conductivity ImG uu (k = 0, uS) in the low- 
est order (second order) the coupling does not contribute. 
This shows that when P c ^ the phonons along the same 
direction do not couple to the magnons. These phonons 
make no contribution to the electromagnons that may be 
observed in an optical conductivity measurement. How- 
ever, the higher order terms are still coupled. A similar 
type of result can be obtained with the terms involv- 
ing SS2- This verifies the first part of our selection rule 
by showing that when an electric polarization is present 
along the c-axis only the phonons polarized along the 
a-axis can couple to the magnons to generate the elec- 
tromagnons. 

The coupling between the phonons polarized along the 
a-axis and the magnons can be derived as follows. For 
u a , from Eq. [38] we have 

u a (5Si x dyS 2 Q - S 2 o x d y SSi) b 
= S 2 u a (-q cos(qy)m^ + s\a{qy)d y rn 2 1 " > ) (50) 

Therefore the total effective Hamiltonian, Hj , describ- 
ing the dynamics of the phonon and the magnon can be 
written as 

PP 2 



p i 

b Hp - q))[-qm+(p) - pm+(p)]}(5i) 



where h is the Planck's constant, p ~ J (the effective 
magnetic coupling strength in the magnetic chain), q is 
the incommensurate wavevector, wo is the bare phonon 
energy, b p is phonon annihilation operator, and m + = 

A qualitative understanding of the electromagnon fre- 
quency for the L1CU2O2 compound can be obtained by 
using a spinwave analysis. The analysis was done for the 
case when the optical phonon frequency is much larger 
than the magnon frequency of interest. This implies that 



if we measure the ac-conductivity versus frequency one 
will observe a peak at the frequency of the magnon. We 
should also note that two points justify the neglect of the 
dynamic lattice degrees of freedom, the displacement field 
u, for the L1CU2O2 system. First, the optical phonon 
frequency is much higher than the magnon frequency. 
Second, the magnetoelectric coupling is very small com- 
pared with to the other multiferroics, like the 113 sys- 
tems, making it unnecessary to explicitly include the dy- 
namic degrees of freedom in the dielectric displacement 
(23|. We now perform a standard spinwave analysis [23[ 
about the ground state to find the frequency at the de- 
sired wavevector, k = Q. 

The model Hamiltonian for L1CU2O2 proposed by [32j 

is 



1,3 



-J^Si^j ■ Si+4j + J±Sij ■ Si.j+i) + DSj_. 



(52) 



Here we have included an easy plane anisotropy which 
generally exists due to the anisotropy of the lattice and 
this term favors the easy plane of spin alignment in the 
ground state, say, the bc-plane in zero field, but the ac- 
plane in a strong magnetic field along the b-axis. We 
solve for the spin wave dispersion. We employ the ro- 
tating frame of reference coordinate system to write the 
spin vector at any site relative to the other rotating one 
as 

S id = Sl d e x + [Si- cos(Q • Ri,j) + sin(Q • Rij)]e y 



+[-S?i sin(g • Ri,j) + cos(Q • R 



(53) 



In the above equation the magnetic propagation vector 
is given by Q = ^,0) where £ = 0.174 is the spiral 
modulation along the chain. Inserting the above form 
for the spin vector we recast the Hamiltonian in terms 
of S^,S n and components. We then compute the 
effective magnetic field components using the formula 
H e ff = ~^§- We then solve for the equation of 
motion for the spin components using the effective mag- 
netic field equations and spin components in terms of the 
Holstein-Primakoff formalism listed below 



(54) 
(55) 

(56) 
(57) 



•5*™ ™ \l ~ {Q"tn,n &rri n) 



^m,n = ~*y ~^( a m,n ~ a m ,ri) 

The equations of motion are given by 

" /J m,n — °m,n 11 ef f °m,n 11 ef f 



l eff 



where in the above equation we set S% j = S since it is 
the direction in which the spin average points. Using 
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the equations for the effective magnetic field, the spin 
components, and finally Fourier transforming we obtain 
the dispersion as 

- ^(q) = SHq)(J(q) J(Q) + D), (59) 

where we define J(q) = J\ cos(q6&) + J 2 cos(2qtb) + 
J4 cos(4qb6) + J± cos(q a a) where q a and qb refer to the 
a and b components of the wavevector. From the equa- 
tions above it is obvious that at q — Q we have 

MQ) - 2S^( Ji2 Q ) + J{0) --,J(Q))D. (60) 

The electromagnon dispersion indicates that the fre- 
quency is proportional to the square root of the easy 
plane anisotropy, D. In zero field, the easy plane is the 
6c-plane. With an applied magnetic field along the a- 
axis, this field will effectively increase the anisotropy in 
the form of D ~ D + H%/2p [33], where p is the spin 
stiffness. This should lead to a hardening of the elec- 
tromagnon frequency (a right shift of the peak in the 
ac-conductivity measurement). If the magnetic field is 
applied along the o-axis, it effectively diminshes the easy 
plane anistropy in the form of D ~ D — H b 2 /2p. There- 
fore, one must observe the softening along with the in- 
crease of the magnetic field and at the point where the 
frequency becomes zero, the magnon mode becomes un- 
stable and the spin-flop transition happens 

The above mechanism explains the observed phase 
transition at a certain magnetic field along the &-axis 
as the destabilization of one electromagnon mode. In 
the high field phase, the spins are lying in the ac- 
plane, making the ac-plane an easy plane described by 
an anisotropy term. In this phase, further increasing 
the field along the &-axis will increase the easy plane 
anisotropy, whereas applying a field along the a-axis will 
decrease the anisotropy. Therefore, we predict that the 
electromagnon hardens with an increase of magnetic field 
along the 6-axis, and softens with an increase of magnetic 
field along the a-axis. This is opposite to what we should 
see in the low field phase. 



VI. SUMMARY AND CONCLUSIONS 

In this paper we analyze the possible types of mag- 
netoelectric coupling in the recently studied multiferroic 
compound LiC^Oi- Based on a group theoretical anal- 
ysis we construct a multi-order parameter phenomeno- 
logical model for the double chain zig-zag structure. We 
show that a coupling involving the inter-chain magnetic 
structure and ferroelectricity is crucial in understanding 
the results of Park et.al. [l(. This constructed model 
for the multiferroic L1C1L2O2 compound can explain the 
polarization flop from the c-axis to the a-axis with the 
applied magnetic field along the &-axis. The model also 
makes specific selection rule predictions about the hybrid 
phonon and magnon excitations called electromagnons. 
We also predict that the electromagnon peaks measured 
in an ac-conductivity measurement are field dependent 
and behave in opposite ways in the P/ / a phase and P/ / c 
phase. However, since the value of the polarization in this 
material is rather weak it will require a very high reso- 
lution spectroscopy measurement to observe the electro- 
magnons in the actual system. 

The results in this paper suggest that the effective the- 
ory of magnetoelectric coupling is richer in a system de- 
scribed by multi-order parameters where several magne- 
toelectric coupling terms can be constructed, than that in 
a system described by a single order parameter. Similar 
physics has also been demonstrated for the RMu^Of, [37| 
multiferroic system and it should be thoroughly investi- 
gated for the other materials. 

The model we propose in this paper could be oversim- 
plified. However, at present only a limited set of exper- 
imental results are available for this compound. We be- 
lieve this phcnomcnological model is a first step towards 
understanding the unique and novel magnetoelectric cou- 
pling observed in the LiCu^O-i compound. It is open 
to future experiments to determine the relevance of the 
other magnetoelectric coupling terms which we derived 
in this paper based on a group theoretical calculation. 
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